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Abstract We assess the statistical spatial distribution of plasma in Mercury's magnetosphere from 
observations of magnetic pressure deficits and plasma characteristics by the MErcury Surface, Space 
ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft. The statistical distributions of proton 
flux and pressure were derived from 10 months of Fast Imaging Plasma Spectrometer (FIPS) observations 
obtained during the orbital phase of the MESSENGER mission. The Magnetometer-derived pressure 
distributions compare favorably with those deduced from the FIPS observations at locations where 
depressions in the magnetic field associated with the presence of enhanced plasma pressures are 
discernible in the Magnetometer data. The magnitudes of the magnetic pressure deficit and the plasma 
pressure agree on average, although the two measures of plasma pressure may deviate for individual 
events by as much as a factor of ~3. The FIPS distributions provide better statistics in regions where the 
plasma is more tenuous and reveal an enhanced plasma population near the magnetopause flanks resulting 
from direct entry of magnetosheath plasma into the low-latitude boundary layer of the magnetosphere. The 
plasma observations also exhibit a pronounced north-south asymmetry on the nightside, with markedly lower 
fluxes at low altitudes in the northern hemisphere than at higher altitudes in the south on the same field line. 
This asymmetry is consistent with particle loss to the southern hemisphere surface during bounce motion in 
Mercury's offset dipole magnetic field. 


1. Introduction 

Magnetometer [Anderson et at., 2007] observations by the MErcury Surface, Space ENvironment, GEochemistry, 
and Ranging (MESSENGER) spacecraft [Solomon et at., 2001] have proven to be one way to characterize the 
plasma pressure distribution in Mercury's magnetosphere. The diamagnetic character of plasmas requires a 
correspondence between plasma pressure and a decrease in magnetic pressure, p B = B 2 /2/j 0 , where B is the 
magnetic field magnitude and // 0 is the magnetic permeability of free space. A localized plasma pressure is 
evident by a depression of the magnetic field. Diamagnetic depressions are frequently observed in the 
MESSENGER Magnetometer (MAG) data within Mercury's magnetosphere and have been used to determine 
the distribution of enhanced magnetospheric plasma pressures [Korth et at., 2011, 2012]. In this paper, we 
extend those earlier analyses with plasma observations. 

The spatial distribution of plasma pressure is important for understanding the interactions among the solar 
wind, magnetosphere, and planetary surface. Plasma enhancements near the magnetopause indicate where 
solar wind plasma enters the magnetosphere via magnetic reconnection [Parker, 1957; Petschek, 1964], 
diffusion [Eastman etal., 1985], and turbulent mixing [Axford, 1964] processes. Magnetic reconnection occurs 
both at the dayside magnetopause, where interplanetary and planetary magnetic field lines merge and are 
subsequently swept anti-sunward by the solar wind, and in the magnetotail, where open field lines connected 
to the northern and southern magnetic poles reconnect pairwise and cycle back toward the dayside [Dungey, 
1961]. In the process, solar wind plasma on interplanetary magnetic field lines is introduced to field lines 
rooted in the planet, via both convection and accelerated flows from dayside reconnection. Because 
Mercury's internal magnetic field is weak — the dipole moment is 1 90 nT R M 3 [Anderson et at., 201 2] and the 
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subsolar magnetopause stands off the solar wind at a distance of, on average, only 0.45 R M from the surface 
[Winslow et al., 2013], where R M is Mercury's radius or 2440 km — it has been proposed that reconnection at 
the dayside magnetopause may erode the subsolar magnetosphere and allow direct access of solar wind 
particles to the planetary surface [ Slavin et al., 2007]. Although MESSENGER has not observed the cusp at 
the subsolar magnetosphere [Winslow et a!., 2012], corresponding to a fully eroded magnetosphere, there 
is plentiful evidence for reconnection through frequent loading and unloading of the magnetic tail on 
~2 min timescales, commensurate with that for the Dungey cycle [Slavin et al., 2010a]. The occurrence and 
rate of reconnection can be inferred from observations of flux transfer events (FTEs), which have been 
detected at the equatorial magnetopause under both northward and southward interplanetary magnetic field 
during the first and second Mercury flybys [Slavin eta!., 2010b], in the magnetotail lobes [Slavin eta!., 2012], and 
in the dayside magnetosphere [Imber et al., 201 3]. Statistical analysis by Imber et al. [201 3] showed that 
FTEs alone account for -30% of the total magnetic flux transfer and are thus a major contributor to 
plasma entry into the magnetosphere. 

In addition to reconnection, plasma can enter the magnetosphere through turbulent mixing, and 
Kelvin-Flelmholtz (KFH) waves have been observed at Mercury's magnetopause [Boardsen et al., 2010; 
Sundberg etal., 2010, 201 1], KFH waves arise at the magnetopause from the velocity shear between the anti- 
sunward solar wind flow and the sunward motion of plasma inside the magnetosphere. It has been 
suggested that the finite ion gyroradius affects the growth rate of the KFH instability [Nagano, 1979]. On the 
dawnside, ion gyration correlates with the vorticity of KFH waves, which promotes the formation of waves 
and results in larger growth rates, and thus larger entry rates, than on the duskside, where ion gyration 
reduces the KFH instability. Flowever, G lassmeier and Espley [2006] predicted the duskside to be more 
susceptible to the KFH instability, and recent particle-in-cell [ Nakamura et al., 2010] and kinetic hybrid 
simulations [Paral and Rankin, 2013] demonstrate that the finite gyroradius of magnetosheath ions and 
the direction of convection electric fields broaden the shear layer at dawn, which lowers the growth rate of 
the KH instability in this region. Consistent with the latter scenario, a preliminary statistical study by 
Sundberg et al. [201 2] identified signatures of KFH waves in the Magnetometer data predominantly at dusk. 
Regardless of preferential occurrence, the KFH instability has been reported at Mercury and is a candidate 
plasma-entry mechanism for Mercury's magnetosphere. 

Within the magnetosphere, the plasma distribution provides insight into the transport of charged particles 
toward the planet and to the surface. Plasmas entering the nightside magnetosphere are transported toward 
the equatorial plane to form the plasma sheet and thence to the dayside by convection [Roederer, 1970]. Prior 
to the availability of extensive in situ observations, the plasma distribution within Mercury's magnetosphere 
had been determined from numerical modeling. Single-particle trajectory calculations by Delcourt et al. 
[2003] revealed a substantial sodium-ion population near the planet on the nightside. These workers also 
found that the gyroradius of these ions exceeds the length scale of magnetic field variations in most of 
Mercury's magnetosphere, so that sodium and other heavy ions follow nonadiabatic trajectories, the 
behavior of which in turn leads to chaotic pitch-angle scattering [e.g., Buchner and Zelenyi, 1 989] and 
precipitation to the surface within a narrow latitude range in each hemisphere. The spatial distribution of 
protons had been inferred with kinetic hybrid simulations [Travnlcek et al., 2007]. These simulations 
suggested that, under idealized conditions, protons can form a closed equatorial belt around the planet. 

Mapping of plasma pressure enhancements inferred from localized depressions in the magnetic field 
magnitude to the magnetic equatorial plane revealed a large-scale plasma population near the planet in the 
portion of the nightside magnetosphere corresponding to the plasma sheet [Korth et al., 2011]. This 
population is present within a half torus that extends from dusk to dawn at radial distances between 1 and 
2.5 R m from the dipole axis and symmetrically about the magnetic equator and displays an inward radial 
pressure gradient [Korth et al., 2012]. The crosstail current sheet is embedded in the plasma sheet, and 
rotations in the magnetic field observed by the MESSENGER Magnetometer indicate that it is ~0.2 R M thick 
[Johnson et al., 2012]. Near the planet, the plasma sheet extends along magnetic field lines to low altitudes at 
mid-latitudes, where charged particles precipitate and contribute to space weathering of the planetary 
surface. On the dayside, plasma associated with the northern magnetospheric cusp has been observed [Korth 
etal., 201 1, 2012; Zurbuchen et al., 201 1; Raines et al., 2013]. Cusp plasma has direct access to the surface, and 
the northern cusp has a mean extent of -1 1° in latitude and 4.5 h in local time [Winslow et al., 201 2]. The area 
of the southern hemisphere cusp, which cannot be observed by MESSENGER because of the eccentricity and 
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high northern periapsis of the spacecraft's orbit, is estimated to be -4 times larger because of the northward 
offset of the internal dipolar field [Winslow et al., 2012], and the flux to the surface is correspondingly higher. 
The northern cusp population has also been shown to be dynamic, exhibiting plasma pressures for anti- 
sunward orientation of the interplanetary magnetic field (IMF) to be 40% higher than for sunward 
orientation [Winslow et al., 2012]. In summary, mapping the plasma distribution inside the magnetosphere 
yields indirect insight into a range of interaction processes among the solar wind, Mercury's magnetic field, 
and the planetary surface. 

Knowledge of the plasma distribution also contributes to determination of the planetary magnetic field. 
The location of the plasma sheet at the magnetic equator provides independent confirmation of the 
substantial, -480 km, northward offset of the internal dipole [Anderson etal., 201 1, 201 2]. The plasma sheet 
thickness can constrain the thickness of the crosstail current sheet [Johnson et al., 2012] in global 
magnetospheric magnetic field models [e.g., Alexeev et al., 2010]. The effects of plasma pressure on the 
locally measured magnetic field must be considered in determinations of higher-order moments of the 
planetary magnetic field. Estimates for the quadrupole term and upper limits on higher-order moments, 
expressed as Gauss coefficients for the axially aligned terms of degrees / = 2, 3, and 4 obtained by Anderson 
et al. [2012], were g 2 o = — 75 nT, |g 30 | < 22 nT, and |g 40 | < 6 nT. At the magnetic equator, where plasma is 
most abundant in the nightside magnetosphere, the magnitude of the quadrupole term is 


where r is the radial distance. Magnetic depressions in excess of 50 nT have been observed at -1 .2 R M radial 
distance from the internal dipole [Korth et al., 201 1 , 201 2], whereas the quadrupole field at this distance is 54 nT, 
comparable to or smaller than the diamagnetic depressions of the magnetic field. The higher-order fields are 
smaller than the quadrupole and fall off more quickly with increasing distance and thus yield signatures at 
MESSENGER altitudes that are lower than many magnetic depressions. Reliable estimates of the higher-order 
structure in Mercury's magnetic field must therefore account for the effects of the plasma on the magnetic 
field measurements. 

The paper is organized as follows. In section 2, we describe a modification to the magnetic field model used 
for tracing the observations to eliminate mapping artifacts in the equatorial distributions. Analysis of plasma 
data and their spatial distributions are described and compared with those inferred from magnetic pressure 
deficits in section 3. The results are discussed in section 4 and summarized in section 5. 

2. Magnetic Field-Line Mapping 

2.1. Magnetospheric Equator 

The statistical plasma distribution in Mercury's magnetosphere was previously inferred from mapping the 
magnetic pressure deficits observed from orbit over an -1 1 month period along the magnetic lines offeree 
to both the magnetic equatorial plane and to low altitudes [Korth et al., 201 2]. This mapping technique 
requires that the plasma pressure be isotropic, which implies that the pressure is constant along the magnetic 
field line [Spence et al., 1989; Goertz and Baumjohann, 1991]. The global magnetic field needed for the 
mapping process was provided by the Alexeev et al. [2010] magnetic field model, which includes the 
planetary field represented by a spin-axis-aligned, southward-directed dipole of moment 190 nT R M 3 and a 
479 km offset northward along the spin axis [Anderson et al., 2011, 2012; Johnson et al., 2012] and external 
contributions to the magnetic field generated by magnetopause and magnetotail currents. The best-fit 
parameters for the Alexeev et al. [2010] model, which minimize the residuals with respect to the magnetic 
field observations, were determined by Johnson et al. [2012]: the magnetopause subsolar stand-off distance 
from the internal dipole is R ss = 1.45 R M ; the half-thickness of the cross-tail current sheet is 0.09 R M ; the 
distance from the planet surface to the inner edge of the cross-tail current sheet is R 2 = 1.41 /? M ; the 
magnetotail field intensity is 6 T = 137 nT; and the paraboloid flaring parameter is y= 1. Unfortunately, this 
parameterization yields unphysical magnetic islands (O-type magnetic field lines, or O-lines, which close on 
themselves) in the model at the dawn and dusk flanks of the magnetosphere near the planet, attributable to 
an overly sharp inner edge of the tail current sheet. These portions appear as white patches in regions of the 
equatorial plane not accessible by field lines traversed by the MESSENGER spacecraft [cf. Korth et al., 2012, 
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Figures 9-11]. Eliminating O-lines in the 
model required an eightfold increase in the 
half-thickness of the current sheet, from 
0.09 R m to 0.75 R M , and this modified model 
is used to compute a revised equatorial 
map. To quantify the corresponding change 
in the mapping, Figure 1 compares the 
magnetic field lines of the best-fit model 
(solid black lines) and the modified model 
with increased current-sheet thickness (dashed 
red lines) in the noon-midnight meridian of the 
Mercury solar orbital (MSO) coordinate system, 
where +X is toward the Sun, +Y is duskward, +Z 
is northward, and the origin coincides with the 
planet center. Both parameterizations yield 
similar field line geometry except near the 
nightside equatorial plane, where the closed 
field lines of the modified model map 
somewhat closer to the planet. 

An estimate of the mapping uncertainty was 
obtained by mapping magnetic field lines 
originating at the planetary surface to the 
magnetic equatorial plane and comparing 
the equatorial crossing points between the 
two models. Figure 2 shows the dependence 
of the equatorial mapping distance on the 
Mercury solar magnetospheric (MSM) 
latitude of field lines originating on the 

planetary surface at local midnight for the best-fit (solid line) and modified (dashed line) Alexeev et at. 
[2010] models. MSM coordinates are centered on the planetary dipole and correspond to MSO 
coordinates displaced 479 km to the north [e.g., Korth et al., 2011]. At low latitudes, the field lines of 
surface locations that increase in latitude map to increasingly distant locations in the equatorial plane. 

Toward mid-latitudes, the equatorial 
mapping point approaches asymptotically 
the latitude of the boundary between open 
and closed field lines. The mapping 
uncertainty associated with the modification 
of the Alexeev et al. [2010] model can be 
estimated from the difference in the open- 
closed boundary between the two models 
and is ~7° in latitude. Similarly, the 
equatorial mapping uncertainty for surface 
locations at a given latitude can be 
estimated from the vertical difference 
between the two curves in Figure 2. This 
uncertainty increases with latitude and 
amounts from <0.5 R M for latitudes <30° to 
~3 R m at the open-closed boundary of the 
best fit model. For altitudes above the 
surface, the equatorial mapping uncertainty 
becomes smaller, so that the estimates 
inferred from Figure 2 are upper limits. 
Furthermore, because the tail current sheet 


-2 0 
X MSO [R m ] 

Figure 1. Comparison of the best-fit (solid black lines) and 
modified (dashed red lines) model [Alexeev et al., 2010] magnetic 
field lines in the MSO X-Z plane at 7=0. The half-thickness of the 
crosstail current sheet for these models is 0.09 and 0.75 ft M , 
respectively. The circle outlines Mercury's surface, and the Sun is to 
the right. 



MSM Surface Latitude [°] 

Figure 2. Dependence of the equatorial mapping distance on the 
MSM latitude of field lines originating on the planetary surface at 
local midnight for the best-fit (solid line) and modified (dashed line) 
Alexeev etal. [2010] models, labeled by the half-thickness of the cross- 
tail current sheet. The asymptotic vertical lines indicate the latitudes 
of the boundary between open and closed magnetic field lines. 
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is restricted to the nightside 
magnetosphere, the latitude difference 
between the open-closed boundaries of 
the best-fit and modified models, and 
thus the equatorial mapping uncertainty, 
decreases toward the dayside and 
vanishes near local noon. 

With the modified Alexeev et al. [2010] 
magnetic field model, we then remapped 
the equatorial magnetic pressure 
distribution [cf. Korth et al., 2012, Figure 9]. 
The revised distribution is shown in 
Figure 3 for a larger section of the magnetic 
equatorial plane than presented in the 
previous study. To account for seasonal 
variations in the solar wind environment 
resulting from the large eccentricity of 
Mercury's orbit about the Sun, all pressures 
were normalized to Mercury's mean 
heliocentric distance, 0.39 All, with the 
scale factor (rj 0.39) 2 , where r s is Mercury's 
heliocentric distance in AU at the time of 
observation [ Korth et al., 2012]. The voids at 
the dawn and dusk flanks are no longer 

present, whereas the dominant features in the equatorial pressure distribution, including the nightside 
enhancement near the planet and the inward radial pressure gradient, are qualitatively preserved. At larger 
distances, especially on the duskside, the distribution becomes patchy where bins are not threaded by field 
lines on which magnetic depressions were observed. The revised distribution in Figure 3 is used below for 
comparison with the plasma observations. 


-2 -1 o 
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Figure 3. Distribution of the mean magnetic pressure deficit mapped 
to Mercury's magnetic equatorial plane and normalized to a heliocentric 
distance of 0.39 AU. Observations were mapped with the modified 
Alexeev et al. [2010] model and a half-thickness of the cross-tail current 
sheet of 0.75 R^. The circle denotes the planet, the Sun is to the right, 
and the magnetopause of the magnetic field model is represented by 
the solid black line. 


2.2. Invariant Latitude 



X MSM [RJ 


Figure 4. Illustration of the mapping procedure to identify the invar- 
iant latitude,/]. Shown in the MSMX-Z plane at 7=0 are the planetary 
surface and the 1 Rm sphere in MSM coordinates (black and grey 
circles, respectively), magnetic field lines (grey), and the MESSENGER 
trajectory in this noon-midnight plane (red, with arrows indicating 
the direction of spacecraft motion). 


In addition to the equatorial distribution, 
the inferred plasma pressures have been 
mapped to low altitudes. To preserve the 
geometry of the field lines on which 
observations were made, the observations 
were mapped onto a grid of invariant 
latitude and local time. The invariant 
latitude, X, is the latitude at which the field 
line threading the observing location maps 
onto a sphere of radius 1 R M centered on the 
offset dipole moment, as illustrated 
in Figure 4. The invariant latitude was 
obtained by tracing the magnetic field line 
threading the spacecraft location to the 1 fi M 
MSM sphere in the respective hemisphere 
with the modified Alexeev et al. [2010] 
magnetic field model. Organizing the 
observations by invariant latitude rather 
than planetocentric latitude facilitates 
critical examination of asymmetries of the 
low-altitude plasma distribution with 
respect to the magnetic equator. 
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Figure 5. Distribution of the mean magnetic pressure deficits with 
respect to invariant latitude and local time observed on the (a) 
descending and (b) ascending legs of each orbit. The observations 
were normalized to a heliocentric distance of 0.39 AU. The black lines 
confine the region of closed model field lines near the equator not 
sampled by MESSENGER. 


The low-altitude distribution of the 
magnetic pressure deficit as a function of 
invariant latitude and local time is shown in 
Figure 5, in which observations on the 
descending (Figure 5a) and ascending 
(Figure 5b) orbit segments are shown 
separately. The equatorial gap in X 
(delimited by black lines in Figure 5) 
corresponds to field lines that are not 
intersected by the MESSENGER orbit. Both 
panels of the figure show a persistent 
plasma population associated with the 
high-latitude dayside cusp region. 
Observations on the descending orbit leg 
are obtained at lower altitudes than for the 
ascending orbit leg and, hence, extend to 
lower |A| than those on the ascending leg. 
The distributions for the ascending and 
descending leg differ most on the 
nightside. On the descending leg of the 
orbit, there is a prominent north-south 
asymmetry in the pressure, such that 
the latitude range of the pressure 
observations is about three times wider in 
the south than the north. This asymmetry 
is not found on the ascending leg, for 
which the latitude range over which the 
pressure deficits are observed is about 
equal in both hemispheres. 


3. Plasma Data Analysis 

Plasma observations are obtained by MESSENGER'S Fast Imaging Plasma Spectrometer (FIPS), one of two 
sensors on the Energetic Particle and Plasma Spectrometer (EPPS) instrument [Andrews et al., 2007]. FIPS 
combines a time-of-flight detector with an electrostatic analyzer to measure the mass per charge ( m/q ) and 
energy per charge (E/q) ratios of incident ions. The instrument has a simultaneous 1 .4 ji sr field of view, which 
is reduced to 1.1 5n sr by obstructions from the spacecraft, and registers incident ions with an angular 
resolution of -15° [Andrews et al, 2007]. We used measurements of the proton flux as function of E/q to derive 
estimates of the proton thermal pressure, and we employed fits to fill in unobserved portions of phase space 
[cf. Gershman et al., 2013, and references therein] to derive average spatial distributions for comparison with 
the magnetic pressure deficits. Protons were detected by FIPS in 50 logarithmically spaced steps in E/q in the 
energy range from 1 00 eV/e to 1 3.3 keV/e with an integration time of 50 ms at each step. The observations 
included in the statistical analysis were obtained from 1 1 April 201 1 to 12 February 2012, which differs from 
that studied by Korth etal. [2012] (23 March 201 1 to 12 February 2012) by about 2 weeks, owing to differences 
in the completion date of commissioning of the MAG and EPPS instruments for the Mercury orbital phase of 
the MESSENGER mission. We restricted the data set to observations obtained within the magnetosphere, 
which is identified from the MAG data as the interval between the innermost inbound and outbound 
magnetopause crossings [Winslow et al., 2013]. The same normalization for heliocentric distance used for the 
magnetic pressure deficit analysis was applied to the proton counts and thermal pressure estimates. 

3.1. Proton Flux Distributions 

The equatorial distribution of the proton flux was calculated by selecting observations on closed model 
magnetic field lines, mapping them to the magnetic equator, and averaging them in bins of dimensions 
0.1 25 x 0.1 25 R m . The resulting distribution in the MSM X-Y plane at Z= 0 R M is shown in Figure 6. The 
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enhanced plasma population observed in a 
toroidal section extending near the planet 
nightside from dusk to dawn is similar to 
that inferred from the magnetic pressure 
distribution (cf. Figure 3). The low flux levels 
on the dayside confirm that the plasma 
population is tenuous in this region, 
supporting the hypothesis that protons in 
the energy range measured by FIPS 
typically do not complete their drift orbits 
around the planet as they do at Earth [Ness, 
1979; Russell et a!., 1988; Korth et at., 2012]. 
Note that the increase in plasma pressure 
at the dayside magnetopause inferred from 
the Magnetometer data was shown to be 
due to cusp pressure observations having 
erroneously been identified as located on 
closed field lines in the model. In contrast 
to the plasma distribution inferred from 
Magnetometer data, Figure 6 also shows 
increases in the proton flux levels at 
magnetopause flanks, which are stronger 
at dawn than at dusk. These enhancements are less pronounced in the magnetic pressure distributions 
because the event count is smaller. 

To illustrate that the features evident in the statistical analysis correspond to structures observed on 
individual orbits, FIPS proton and Magnetometer data are shown for four representative magnetospheric 
transits in Figure 7. The innermost inbound and outbound magnetopause boundaries are indicated by the 
vertical dashed lines. In Figure 7c, the proton spectrogram changes near 2100 UTC, when the FIPS instrument 
transitioned to a scan mode with lower time resolution and lower minimum energy. Figures 7a and 7b show 
data from orbits with spacecraft motion from dawn to dusk and dusk to dawn, respectively. During both 
orbits, enhanced proton fluxes were evident at the dawnside and duskside magnetopause flanks adjacent to 
the magnetopause, which is unambiguously identified by a rotation in the magnetic field azimuth angle and 
a change in the amplitude of magnetic fluctuations. In addition, the magnetopause distinctly separated the 
more energetic protons observed at the magnetosphere flanks from the colder magnetosheath population. 
For both cases in Figures 7a and 7b, the count rates at dawn were substantially higher than those at dusk, 
consistent with the statistical dawn-dusk asymmetry and indicating that the asymmetry is at least partially 
free of bias in the observations with respect to the Mercury's orbit phase. The flank populations were not 
previously identified in the Magnetometer observations. Comparison of depressions in the magnetic field 
magnitude associated with the flank populations (Figure 7, top panels) with those observed in the plasma 
sheet and cusp [Korth et al., 2011, 201 2] shows that the former are small compared with the latter and occur 
over extended periods of time so they are difficult to identify reliably in the Magnetometer data. 

Observations during magnetosphere transits with the spacecraft moving in midnight-to-noon and noon- 
to-midnight directions are shown in Figures 7c and 7d, respectively. The key features observed in these orbits 
are as follows. The flux levels are low on the dayside equatorward of the cusp, 2007 to 2015 UTC on orbit 205 
and 2214 to 2229 UTC on orbit 500, demonstrating that few particles complete their drift around the planet 
dayside. The fluxes are also low at middle to high latitudes, 2021 to 2031 UTC (northern hemisphere) and 
2052 to 21 29 UTC (southern hemisphere) on orbit 205 and 2030 to 21 07 UTC (southern hemisphere) and 2202 
to 2207 UTC (northern hemisphere) on orbit 500, corresponding to the region that is magnetically connected 
to the center-tail minimum observed between the flank populations in Figure 6. As the spacecraft crossed the 
magnetic equator near midnight, there was a dropout in the proton fluxes north of the equator, from 2033 to 
2038 on orbit 205 and from 2140 to 2145 on orbit 500. 

The low-altitude mapped distribution of the proton flux is obtained in the same manner as that for the 
magnetic pressure deficit by binning the observations with respect to magnetic local time and invariant 
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Figure 6. Distribution of the mean proton flux mapped to Mercury's 
magnetic equatorial plane and normalized to a heliocentric distance 
of 0.39 AU. The circle denotes the planet, the Sun is to the right, and 
the magnetopause of the magnetic field model is represented by the 
solid black line. 


KORTH ETAL. 


©2014. American Geophysical Union. All Rights Reserved. 


2923 


^AGU Journal of Geophysical Research: Space Physics 


10.1 002/201 3JA01 9567 



Figure 7. Magnetic field and plasma observations during magnetosphere transits for the orbits indicated at the top of each panel. From top to bottom, each panel 
shows the magnetic field magnitude; the magnetic field polar (black) and azimuth (red) angles; the proton flux spectrogram color-coded by count rate and 
over-plotted in grey with the total count rate; and the spacecraft MSM latitude (black) and local time (red). The innermost magnetopause boundaries are 
indicated by vertical dashed lines. The panels show sample orbits with the ascending orbit leg in the (a) dawn-to-dusk, (b) dusk-to-dawn, (c) noon-to-midnight, 
and (d) midnight-to-noon directions. 


latitude as described in section 2.2. The distributions for the descending and ascending orbit leg observations 
are shown in Figures 8a and 8b, respectively. Each panel shows the mean proton flux color-coded by 
magnitude in bins of dimensions 0.5 h in magnetic local time and 1° in invariant latitude. The heavy black 
lines in both panels show the mean X closest to X = 0° sampled on the descending orbit leg. The equatorial 
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edge of the distributions corresponds to 
the minimum |/| sampled at that local 
time, that is, the magnetic equator. The 
white and red traces in the figure are 
discussed in section 4. 

The highest proton fluxes are observed on 
the nightside at middle southern latitudes 
and in the northern cusp on the dayside. 
On the ascending orbit leg, the highest 
fluxes are in the northern cusp. In the 
descending leg observations, there is a 
pronounced north-south asymmetry in 
the nightside flux distribution. The most 
intense nightside fluxes are observed in a 
substantially wider latitude range in the 
southern hemisphere than in the north. 
Whereas the southern hemisphere mid- 
latitude flux is uniformly distributed in 
latitude within -20° in X of the equator, in 
the northern hemisphere, there is a band 
of enhanced fluxes within ~5° of the 
minimum X sampled, near the equator, 
and a weaker enhancement in the fluxes 
1 5° to 20° in X northward. As discussed 
below, the north-south asymmetry in the 
nightside descending leg observations is 
attributed to the northward offset of the 
planetary magnetic dipole. 

Both the ascending and descending 
observations show that the fluxes at dawn 
are higher than at dusk, consistent with the 
equatorial mapping (cf. Figure 6). The 
dawn-dusk asymmetry is most evident in the southern hemisphere descending observations, but it is present 
also in the ascending orbit data in both the northern and southern hemispheres. 

Finally, some flux observations map to the dayside polar region but are in fact observed at spacecraft 
locations tailward of the dawn-dusk terminator. These fluxes are observed in close proximity to the 
magnetopause on open magnetospheric field lines, which are rooted in the planet at high latitudes on the 
dayside and extend far into the magnetotail, as shown in Figure 1. The mapping along field lines close to 
the magnetopause is very sensitive to the geometry of the magnetic field and small deviations between 
the model and the natural system and result in large uncertainties, so the results in this region are not regarded 
as indicative of the low-altitude fluxes. Therefore, the southern high-latitude observations are not considered 
further in this study. 

3.2. Proton Pressure Distributions 

The FIPS instrument observes only partial distribution functions, so plasma moments, including the 
thermal pressure, cannot be directly computed from the observations under all conditions. However, for 
plasmas with subsonic flow velocities and modest temperature anisotropies, conditions that are typically 
encountered throughout Mercury's magnetosphere and, particularly, in the nightside plasma sheet, 
density, temperature, and pressure estimates can be recovered from E/q distributions only [Cershman eta!., 
2013]. The above conditions are not as well satisfied in the cusp, where precipitation to the planetary 
surface can result in strong directional anisotropies, but FIPS observations there should nonetheless 
provide a reasonable first-order estimate of plasma pressure. The plasma pressures were derived from 


a) 90 



0.12 


Figure 8. Distribution of the mean proton flux with respect to invar- 
iant latitude and local time observed on the (a) descending and (b) 
ascending legs of each orbit. The observations were normalized to a 
heliocentric distance of 0.39 AU. The black lines confine the region of 
closed model field lines near the equator not sampled by MESSENGER. 
The boundaries between open and closed field lines in the best-fit and 
modified Alexeev etal. [2010] magnetic field model are represented by 
white dashed and solid lines, respectively. The red line in Figure 8a 
depicts the loss boundary. 
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Figure 9. Distribution of the mean proton pressure with respect to invariant latitude and local time observed on the 
(a) descending and (b) ascending nodes of the orbit. The observations were normalized to a heliocentric distance of 
0.39 AU. The black lines confine the region of closed model field lines near the equator not sampled by MESSENGER. 


proton E/q distributions by application of the technique of Gershman et al. [2013] to data acquired over 
1 min intervals within Mercury's magnetosphere during which more than 50 ion counts were recorded, 
summed over all energy channels. The pressure estimates are for protons only and do not account for 
contributions from electrons or heavy ions. 

Statistical distributions of the proton pressure with respect to magnetic local time and invariant latitude were 
obtained in the same manner as for the distributions of magnetic pressure deficit (Figure 5) and the proton 
flux (Figure 8). Because of the ~1 min integration time for each pressure inversion, a larger bin size of 1 h by 
2.5° was used to average the results. The average plasma pressure distributions for the descending and 
ascending legs of the MESSENGER orbits are shown in Figure 9. The distribution of the proton pressure is 
similar to that of the magnetic pressure deficit and the proton flux. The highest pressures are primarily in the 
plasma sheet and the cusp. Furthermore, the nightside plasma sheet pressure observed on the descending 
leg of each orbit is markedly higher in the south than in the north for invariant latitudes within -20° of the 
magnetic equator. There is also a faint double-banded structure in the pressure on the descending leg at 
night in the northern hemisphere, with enhanced pressures within a few degrees of the equator and again 
-20° farther poleward. The pressure distributions for both the descending and ascending orbit segments 
have higher pressures at dawn than at dusk. 

To quantify the comparison between the FIPS proton pressure estimates and the magnetic pressure deficits, 
we paired each plasma pressure estimate with the average of the magnetic pressure deficit derived for the 
same interval. The observations of FIPS-derived proton pressure versus magnetic pressure deficit are shown 
in Figure 10. Figure 10a includes proton pressure observations within the magnetosphere but excluding 
those on the dayside northward of 45°N latitude, the region broadly associated with the northern cusp; the 
observations for the latter region are shown in Figure 10b. The regression coefficients for the sets of data 
points in Figures 10a and 10b are 0.25 and 0.48, respectively, indicating only modest correlation between 
the pressure estimates. From the distributions, it is evident that order of magnitude differences between 
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the Sp B and p H+ values are not 
uncommon and that the dayside high- 
latitude observations exhibit somewhat 
more scatter. The distribution of data 
excluding the cusp peaks for p H+ = Sp B 
and clusters within a factor of 3 of the 
Ph+ = Sp B line. The data in the cusp region 
distribution are biased to p H+ < 6 p B . 

4. Discussion 

Statistical distributions of the proton flux 
and pressure derived from FIPS 
observations in orbit about Mercury 
have been compiled by averaging 
measurements projected along 
magnetic field lines both to the 
magnetic equatorial plane and to low 
altitudes. The FIPS results have also been 
compared with estimates of the total 
plasma pressure from magnetic pressure 
deficits. The FIPS distributions are 
qualitatively similar to those of the 
magnetic pressure deficit in regions 
where magnetic depressions were 
identifiable. In addition, the proton 
pressures derived from the FIPS data and 
the magnetic pressure deficits are in 
agreement on average, although 
substantial differences were found from 
case to case. We first discuss possible 
reasons for the scatter in the 
quantitative comparison of proton 
thermal pressure and magnetic pressure 
deficit results and then proceed to a 
consideration of the key results of the 
plasma distribution in Mercury's magnetosphere and the implications for plasma entry, transport, and 
precipitation loss to the surface. 

Although the distributions of the magnetic pressure deficit and the proton thermal pressure compared well 
qualitatively, direct comparison of the pressures reveals only a modest correlation. The FIPS proton pressure 
estimates agree on average with the magnetic pressure deficits except in the cusp region, where p H+ is 
systematically low. In addition, discrepancies in the pressures of up to an order of magnitude were observed, 
but outside of the cusp region the differences within the magnetosphere were typically within a factor of 3. 
Several factors may contribute to these disparities. First, the magnetic depressions reflect a magnetic 
pressure deficit relative to the background pressure. If the background plasma pressure is not negligible, the 
magnetic pressure deficit underestimates the total plasma pressure. Second, the FIPS pressures are for 
protons only and do not include contributions from electrons or heavy ions. In the Earth's magnetotail plasma 
sheet, the proton temperature is a factor of 8 greater than the electron temperature on average [Baumjohann 
eta!., 1989]. So unless the relative dynamics forming the plasma sheet at Mercury are markedly different from 
those at Earth, ignoring the electron pressure should lead to a modest, 10 to 20%, underestimate of the 
plasma pressure. Compositional observations in Earth's plasma sheet have revealed a relative FHe 2+ abundance 
of 3% [ Lennartsson and Shelley, 1 986; Lennartsson, 1 992; Fuselier et at., 1 999], consistent with a relative 
abundance for these ions of 3-5% in the solar wind [ Feldman et al., 1 978] and solar wind penetration into the 
magnetosphere. The temperature of FHe 2+ typically exceeds that of Fl + by a factor of 4 [Fuselier et al., 1 997]. 
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Log Magnetic Pressure Deficit [nPa] 

Figure 10. Proton thermal pressure versus simultaneously observed 
magnetic pressure deficit for (a) non-cusp and (b) cusp events. The cusp 
is broadly defined as the dayside region northward of 45°N latitude. The 
black lines indicate a slope of unity for reference. 
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Thus, if the solar wind heating and acceleration processes at Mercury are comparable with those at Earth, 
He 2+ may contribute as much as 10-15% to the total pressure. In addition, FIPS observations in Mercury's 
plasma sheet during three orbits recorded in April 201 1 have shown contributions from sodium ions with 
pressure ratios of p Na + /p H + = 4-33% [Zurbuchen et al., 201 1], Thus, pressures from other species could 
contribute up to half of the total pressure, so the FIPS proton pressures could be as much as a factor of 2 
lower than <Sp B for this reason. Finally, we expect that the protons in the cusp are highly anisotropic since 
the fluxes directed upward from the planet should be markedly depleted for pitch angles within the loss 
cone [cf. Winslow et al., 2012]. This direction dependence violates the assumption of isotropy applied in the 
FIPS proton pressure inversions and would lead to an underestimate of the proton pressure [Gershman 
et al., 2013]. It is also possible that other species contribute significantly to the plasma pressure in the cusp 
[cf. Raines et al., 201 1, 2013 ; Zurbuchen et al., 201 1], Thus, we conclude that the systematic difference between 
Sp B and p H+ in the cusp is likely due to the character of the proton distributions, whereas the factor of 2-3 
scatter elsewhere could be due to plasma composition variations so that p H+ underestimates the total 
plasma pressure, background plasma effects so that Sp B underestimates the total pressure, and uncertainties in 
the FIPS pressure inversions. 

Nonetheless, the direct measurement of protons by FIPS provides substantially greater fidelity in the 
measurement of the distribution of plasma in Mercury's magnetosphere than the magnetic pressure 
deficits alone. Both the FIPS results and the pressure deficit analysis show strong enhancements in the 
northern cusp region and a plasma enhancement in a toroidal section extending at 1-2 R M radial distance 
on the planet's nightside from dusk to dawn. In addition to confirming these results, however, the FIPS 
observations reveal features that were either difficult to discern or could not be resolved with the 
Magnetometer data. 

The FIPS proton fluxes exhibit a pronounced dawn pressure enhancement along the magnetopause flank, a 
minimum in the center of the far magnetotail, and another enhancement along the dusk flank but with fluxes 
markedly lower than at the dawn flank (Figure 6). This feature of the plasma distribution in the tail was barely 
perceptible in the magnetic pressure deficit analysis. The open-closed boundary (cf. Figure 8) was computed 
by tracing magnetic field lines from the northern to the southern hemisphere surface and vice versa in 0.1° 
latitude steps for each hour in local time until either the opposite hemisphere surface was reached, 
corresponding to closed field lines, or the field line distance exceeded 10 R M , at which point the field line was 
assumed to be open. Thus, invariant latitudes equatorward of these lines on the nightside map to the 
magnetotail, as shown in Figures 3 and 6. The dawn-dusk asymmetry is also prevalent in the low-altitude 
distributions for the southern descending leg and for the ascending leg in both the north and south 
(Figure 8). Recall that the proton fluxes and pressures have been scaled by 1/r s 2 to correct for variations in 
average solar wind background density with Mercury heliocentric distance. The fact that this dawn-dusk 
asymmetry was observed for individual cases at Mercury aphelion (Figure 7c) and perihelion (Figure 7d) 
indicates that the statistical result reflects an actual feature of the plasma distribution at Mercury and is not a 
residual artifact in the seasonal coverage provided by the MESSENGER orbit. 

Plasma enhancements at magnetopause flanks have also been observed in the low-latitude boundary layer 
of the terrestrial magnetosphere [Wing and Newell, 1 998], where particles from the magnetosheath can enter 
the magnetosphere, e.g., through turbulent mixing of the solar wind and the outer fringes of the 
magnetosphere resulting from the KFH instability [Axford, 1964]. Diffusive particle entry into the 
magnetosphere has been observed at Earth [Terasawa eta!., 1997; Fairfield et al., 2000], and the considerable 
mass transport into the magnetosphere has been linked to magnetic reconnection in KFH vortices [Offo and 
Fairfield, 2000], The observed dawn-dusk asymmetry in the proton fluxes is consistent with the hypothesis by 
Nagano [1979] that ion gyration promotes the formation of KFH waves on the dawnside, leading to larger 
growth rates and thus larger entry rates at dawn than at dusk, where ion gyration reduces the KFH instability. 
This hypothesis is, however, inconsistent with observations by MESSENGER [Boardsen et al., 2010; Sundberg 
et al., 2010, 201 1] and particle simulations [Nakamura et al., 2010; Paral and Rankin, 2013], which show KFH 
waves primarily at the duskside magnetopause, where plasma enhancements should thus be stronger. The 
above considerations are predicated on the assumption that the seed plasma population in the 
magnetosheath is relatively uniform, so KFH mixing introduces similar amounts of plasma at both flanks. 
FHybrid simulations [Travnlcek et al., 201 0], however, have shown pronounced dawn-dusk asymmetries in the 
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magnetosheath plasma distribution, which depend on IMF orientation [Herclk et al., 2013]. Thus, whereas the 
KFH instability remains a possible means for diffusive plasma transport into the magnetosphere, additional 
analysis accounting for the magnetosheath source plasma is needed to determine definitively the 
importance of this mechanism over other processes. 

The equatorial distribution of the proton flux (Figure 6) also shows that plasma is depleted in the dayside 
magnetosphere. This result is in contrast to that from the magnetic pressure deficit analysis (Figure 3), which 
did not resolve this feature. Indeed, as noted in the previously published analysis of the magnetic pressure 
deficits [Korth et al., 2012], the large magnetic pressure deficits on the dayside are attributed to magnetic 
depressions within the magnetosheath that are erroneously mapped to the equatorial plane with the 
average magnetopause and magnetic field model. Note that for consistency with previous magnetic field 
analysis, the magnetic pressure deficit statistics were calculated from data intervals constrained by the model 
magnetopause because the magnetopause crossing list was not available when these analyses were first 
done. The FIPS analysis, however, uses the actual magnetopause crossings so that magnetosheath data are 
excluded by definition. Moreover, the dayside magnetosphere is sampled by MESSENGER near the planet at 
fairly low altitudes, where the gradient in the planetary field along the MESSENGER trajectory is high. 
Discerning deficits in the magnetic field against a steep gradient in the observed field is not straightforward, 
but the direct detection of protons by FIPS alleviates this problem. The FIPS observations are thus critical in 
correcting this limitation of the pressure deficit analysis, and they confirm that the plasma typically does not 
encircle the planet on closed drift paths, consistent with Mercury's small dipole magnetic and corotation 
electric fields [Ness, 1 979; Russell et al., 1 988; Slavin et al., 2007; Korth et al., 2012]. Our statistical results do not 
contradict the findings from previous kinetic hybrid simulations, including the prediction of a closed-belt 
proton population encircling the planet near the equator. This feature was most evident in simulations with 
low and constant solar wind pressure [Travnlcek et al., 2007] or northward IMF conditions [Travnlcek et al., 201 0]. 
Under these low-activity conditions, the magnetopause standoff distance, and the size of the magnetosphere in 
general, is larger than encountered for the average magnetosphere. Mercury's space environment is highly 
dynamic [ Slavin et al., 2007, 2010a, 2010b, 2012; DiBraccio et al., 2013], and solar wind dynamic pressure 
and IMF magnitude and direction, which control magnetospheric configuration and dynamics, vary on 
short timescales. In contrast, solar wind conditions that are stable over prolonged periods of time, as 
assumed for the simulations cited above, are rarely observed. Thus, it is not surprising that the observed 
average distribution, which reflects the full range of dynamic conditions, deviates from the simulation results. 

Finally, structure in the nightside plasma sheet with invariant latitude is unmistakable in the FIPS proton data. 
At lower altitudes, the plasma populations observed by FIPS and inferred from the magnetic field data on the 
descending orbit leg show a pronounced north-south asymmetry (cf. Figures 3a, 8a, and 9a). On the 
descending leg of the orbit, the nighttime northern hemisphere proton fluxes are observed in two latitudinal 
bands, a few degrees wide in invariant latitude separated by a gap of 10° to 20° (cf. Figure 8a). This double- 
banded structure is also recognizable in the distribution of the proton pressure (Figure 9a). The southern 
hemisphere distribution exhibits higher fluxes and has no gap. Observations from the ascending orbit leg do 
not exhibit this asymmetry but do not sample as low in invariant latitude owing to the higher altitudes on the 
ascending leg of the MESSENGER orbit during the time the data were acquired for this study. 

There are several features of the invariant latitude distributions that provide clues to the origin of this north- 
south asymmetry in the observed proton fluxes and pressure. First, we note that the high proton counts and 
pressures in the southern hemisphere of the descending leg data (Figure 8a) occur on closed field lines, that 
is, equatorward of the white lines denoting the open-closed boundary. In the north, the increase in proton 
fluxes with increasing A on the nightside occurs near the open-closed boundary. The midlatitude region 
featuring enhanced proton fluxes is largely contained within the closed field line region, where these 
particles bounce along the field lines between hemispheres and can be at least partially trapped. On the 
other hand, the plasma population on open field lines, which do not support trapping, is tenuous. Second, 
the ascending orbit leg observations (Figure 8b) exhibit comparable proton fluxes near the open-closed 
boundary in both the north and south. Third, the northern hemisphere fluxes on the ascending leg are higher 
than they are on the descending leg even though the orbits on average sample the same field lines near the 
open-closed boundary. In the south, in contrast, the ascending leg fluxes are essentially the same as those on 
the descending leg at the same A and local time. Thus, the asymmetry is prevalent primarily or only in the 
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descending node observations, and there is a 
difference in the flux magnitude on the same 
closed field lines when observed on the 
ascending and descending orbit legs in the 
northern hemisphere. 

The observations therefore suggest that 
something must be distinctive about the region 
of space traversed by the spacecraft in the 
northern hemisphere on the descending orbit 
leg such that the proton fluxes are reduced 
relative to their levels elsewhere on the same 
closed field lines. The red line in Figure 8a in the 
northern hemisphere shows the boundary 
where, on average, the local magnetic field 
magnitude observed in the northern 
hemisphere is equal to the southern hemisphere 
surface magnetic field magnitude as calculated 
from either the best-fit or modified Alexeev et al. 
[2010] models. (Both versions of the model give 
the same result for this calculation since the 
thickness of the cross-tail current has little effect on the model field at the planetary surface). For / northward 
of the red line, the local magnetic field at the MESSENGER spacecraft is greater than the southern hemisphere 
surface field on the same field line. For / southward of the red line, in contrast, the in situ magnetic field is 
lower than the conjugate southern hemisphere surface field. Remarkably, the boundary between in situ 
magnetic fields higher and lower than the conjugate southern hemisphere surface field very closely marks 
the equatorward band of enhanced fluxes in the northern hemisphere descending-leg distribution. Where 
the in situ magnetic field exceeds the conjugate southern hemisphere surface field, the fluxes are depressed. 

This result motivates the interpretation illustrated in Figure 11. Charged particles on these magnetic field 
lines will gyrate about the field line and bounce between mirror points depending on the pitch angle 
sin 2 a min = B eq /B m , where the subscripts on field magnitude denote the magnetic equator and the mirror 
point, respectively [e.g., Baumjohann and Treumann, 1 997]. Particles with pitch angle a min = 90° mirror at the 
equator. On any closed field line, the minimum pitch angle, a min , of particles that do not encounter the 
planetary surface is given by sin 2 a min = B eq /B sur f_ s where B sur f,s is the field magnitude at the southern 
hemisphere surface on the field line. The projection of this surface in the northern hemisphere is shown by 
the grey circle labeled "reflected planetary surface" and is obtained by reflecting the planetary surface about 
the magnetic equator. The portions of field lines in the northern hemisphere that lie below this surface are 
populated only by particles for which the mirror points lie below the surface in the southern hemisphere. 
Thus, the region occupied by trapped plasma lies above this surface and is indicated by orange shading in the 
figure. We therefore interpret the region of depressed fluxes in the northern hemisphere on the descending 
orbit leg as indicating passage of the spacecraft through the void region of closed field lines, below the 
population of stably trapped plasma. Protons in this region are observed either on their first bounce after 
entering the closed field region or as they bounce from the southern hemisphere to the north after being 
pitch-angle scattered into the southern hemisphere loss cone. Thus, one expects much lower fluxes in this 
region and perhaps an increase of flux near the open-closed field line boundary, as observed. The result 
implies that plasma sheet precipitation occurs primarily in the southern hemisphere and that this loss process 
produces a substantial loss cone in the plasma sheet proton distribution, which in turn will promote local ion 
plasma-wave instabilities, which are also observed [Boardsen et al., 2009, 2012]. 

5. Summary 

We have extended the analysis of the plasma distribution in Mercury's magnetosphere by combining direct 
observations of protons by the FIPS sensor with total plasma pressures inferred from observations of 
magnetic pressure deficits by the Magnetometer both on the MESSENGER spacecraft in orbit about Mercury. 


Reflected Planetary 
Surface 



Figure 11 . Formation of a plasma void at low altitudes in the 
northern hemisphere of the planet (filled grey circle) with the 
planetary and magnetic equators marked by dashed and solid 
black lines, respectively. The heavy black lines indicate a set of 
closed field lines originating from a spin-axis-aligned dipole 
offset northward as illustrated. See text for details. 
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The distribution of proton fluxes and pressures compares favorably with regions exhibiting depressions in the 
magnetic field associated with the presence of plasma. On average, the magnitudes of the magnetic pressure 
deficit and of the proton pressure compare well, although the plasma pressure may differ for individual 
events by as much as a factor of ~3. The FIPS proton distributions also provide better statistics in regions 
where the plasma is more tenuous or where magnetic depressions are difficult to identify, such as near the 
dayside equator. Moreover, the proton data reveal features not previously identified, including very low 
fluxes in the dayside closed-field-line region, consistent with the absence of closed drift paths at Mercury, and 
an enhanced plasma population near the magnetopause flanks resulting from direct entry of magnetosheath 
plasma into the low-latitude boundary layer of the magnetosphere. The fluxes of protons at the dawn flank 
are markedly higher than at the dusk flank, suggesting that direct entry is greater at dawn, perhaps due to a 
corresponding asymmetry in Kelvin-Helmholtz wave-driven diffusion. Finally, the plasma observations also 
confirm the north-south asymmetry of the planetary dipole magnetic field, which is evident in the proton 
data via a depression in fluxes on closed field lines in the northern hemisphere when the spacecraft samples 
locations where the local field exceeds the conjugate southern hemisphere surface field. This outcome 
implies that particle loss to the southern hemisphere surface during bounce motion in Mercury's offset dipole 
magnetic field is a prominent loss mechanism for charged particles in the plasma sheet. 
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